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Abstract: A simple method for determining the thermal component of the 
0\ ■ EOS of solids under high pressure is proposed. Applications to the interior of 

r-| . the Earth gives results in agreement with recent geophysical data 

o 

^ ■ 1. Introduction 

Imagine a specimen of any chemical substance. In the case of a thermo- 
rS \ mechanical system, its equation of stae (EOS) will have the general form 

c3 ■ f{p,V,T) = 0, where p,V,T denote,respectively the pressure to which the 

system is subjected, its volume and temperature. 

Proposing a realistic EOS for a given class of physical systems is, generally 
speaking a highly non-trivial problem (see,for instance, Eliezer et al. ,1986; 
Schatzman and Praderie,1990 for details ).One usually starts from an as- 
sumed interparticle potential, determines the form of the Hamiltonian,and 
from it the thermodynamical functions and the EOS. In order to render 
the calculations more tractable, .the EOS is usually determined in form of 
isoterms in the p — V plane. The thermal component is often introduced at a 
later stage (for example, Holian,1986;Eliezer et al. ,1986). 



In laboratory studies, EOS of various systems can be determined exper- 
imentally (examples of recent reviews concerning the subject are Jayara- 
man,1986;Jeanloz,1989;Drickamer, 1990). The situation is much more compli- 
cated in astrophysics,because planetary and stellar structure is unaccessible 
to direct experiments. What is observable, however, are the consequences on 
the surfaces and in the vicinity of these objects of the processes occuring in 
their interiors. Progress in the understanding of the Earth's interior has,for 
example,been provoked by a combination of observation in geophysics and 
seismology, with laboratory studies of the relevant materials performed under 
high pressure and high temperature (Jeanloz, 1990a). 

The purpose of this letter is to propose a simple method for the deter- 
mination of the thermal component of the EOS of solids subdued to high 
pressure. Details of the calculations, and and application to the interior of the 
Earth are presented in the following section, while the third part contains a 
discussion of various factors influencing our results. as well as the possibilities 
of their improvement. Physically,our calculations are based on combination 
known results of solid-state physics with a particular semiclassical theory of 
dense matter proposed by Savic and Kasanin (1962/65;for a recent applica- 
tion, see Celebonovic, 1990b, and references given there). Compared with recent 
work on the subject (such as Renero et a/.,1990;Kumari and Dass,1990a,b), 
the approach proposed in this letter has the advantage of physical and math- 
ematical simplicity,but a disadvantage of being applicable only to solids. 

2. Calculations 



One of the general characteristics of solid bodies is that the atoms (or 
ions) in them,preform small vibrations about their equilibrium positions. Let 
A^ be the number of molecules in the body,and denote by z/ the number of 
atoms per molecule. The total number of atoms in the system is Nu and 
the number of vibrational degrees of freedom is SNu (speaking eaxctly,it is 
3A^z/-6, but 3iVz/>6 ). 



A system of 3iVz/ mutually independent vibrational degrees of freedom is 
equivalent to an ensemble of independent oscillators. The free energy of such 
an ensemble can be expressed as 

F = iVeo + ksTY^Hl - eM-P^^a)) (1) 

a 

where the first term on the right side represents the interaction energy of 
all the atoms in the system in the equilibrium state, and the summation is 
carried over all the normal vibrations, indexed by a . 

At low temperatures,only low frequency terms (i.e., sound waves) with 
ksT = hua have an important contirubution to eq.(l).It can be shown (for 
example Landau and Lifchitz,1959) that the number of vibrations dN per 
interval of frequency du is 

dN _ 3Vuj^ 

'd^~2^^ ^ ^ 

The volume of the system is denoted by V,and u is the mean velocity of 
sound. Changing from summation to integration in eq.(l),one finally obtains 

F = Neo + -4=^ / ln(l - exp{-hujJkBT))du (3) 

Ztt^u'^ Jo 

If we apply the standard thermodynamical identity E = F — T^ it 
follows from eq. (3) that the energy of the system is 

E = Neo + Vn\kBT)yiO{m'') (4) 

The mean velocity of sound can be approximated by the Bohm-Staver 
formula 

Il2 = (1/3) (^Zm/Mvl) (5) 

(Bohm and Staver,1951;Ashcroft and Mermin, 1976), where m and M 
denote, respectively,the masses of electrons and ions in the solid,i)i7' is the 
Fermi velocity and Z is the charge of ions. 



A detailed discussion of the basic ideas of the theory of Savic and Kasanin 
(the SK theory for short) has recently been published (Celebonovic,1989d). 
Within the framework of this theory,the energy per unit volume of a solid is 
given by 

E = 2e^Z{NAp/A)*^^ (6) 

where A is the mass number and A^^ denotes Avogadro's constant. 
Using equations (4)-(6), after some algebra one obtains the following ex- 
pression for the temperature of a solid as a function of its density p 

T = 1.4217 * lO^{p/Ay/^\m/Mf/^Z^/'^ (7) 

The thermal component of the EOS of a solid can be obtained by mul- 
tiplying this result by the density (expressed in suitable units). 

As an astrophysical test,eq.(7) was applied to the model of the Earth 
discussed previously within the SK theory (Savic, 1981 and earlier work) .The 
following results were obtained 



Depth (km) 


0-39 


TABL] 
39-2900 


E I 

2900-4980 


4980-6371 


Pmax(g Cm-3) 


3.0 


6.0 


12.0 


19.74 


Z 


2 


3 


3 


4 


T^UK) 


1300 


2700 


4100 


7000 



where the values of the temperatures have been rounded to the nearest 
±100K. 



3. Discussion 



The distrribution of temperature with depth within the Earth,or any 
other celestial body,is not directly measurable. In order to draw conclusions 
about the thermodynamics of our planet's interior, geophysicists are bound 
to combine examination of rock samples from the outher ~ 200 km of the 
Earth, with high pressure-high temperature experiments (Jeanloz,1990a,b).It 
has thus been shown that, for example,T = (4500 ± 500)7^ at the base of the 
mantle, and that T = (6900 ± 1000)i^ in the center (Jeanloz, 1990a). 

These experiments were performed on materials known to exist in a thin 
layer beneath the Earth's surface, and assumed to represent its composition 
in bulk. One could conjecture that, such an assumption not being directly ver- 
ifiable by experiments, it has strong influence on the results. Some influence 
it certainly has,but it can not be extremely important. Namely,it is possible 
within the SK theory,to determine the bulk chemical composition of an as- 
tronomical object {i.e., the mean value of the mass number of the mixture 
of materials that it is made of);the only input data needed for such a calcu- 
lation are the mass and radius of the object. It turns out that the value of A 
obtained in this way corresponds closely to the value of A for the materials 
used in experiments. 

The calculations discussed in this letter are strongly influenced by a com- 
bination of several factors from the domain of solid state physics. 

The upper limit of integration in eq.(3) is infinite, which simplifies mathe- 
matics,but renders physics unrealistic - there are no infinite frequencies. The 
difficulty could be, at first sight, circumvented by introducing a suitable cut-ff 
frequency,in the spirit of Debye's model. However, one would then encounter 
the following obstacle, which would be the density dependence of the cut-off 
frequency (and the coresponding cut-off temperature). Solving this problem 
demans a knowledge of the elastic constans and inter-particle potentials of 
the material ( a solution for some types of crystal lattices within Debye's 
model is given in de Launay,1953,1954). 



It is clear that attempting to perform such a calculation for the case of 
the interior of the Earth would quickly render the results questionable, due to 
the accumulation of various approximations. Another "solid-state" factor in- 
fluencing our results is the method of calculation of the speed of sound. The 
use of the Bohm-Staver formula amounts to taking into account only the 
main term in the calculation of the band-structure energy of a solid (see Har- 
rison, 1989 for details), which is ,in turn, used in determining the value of the 
velocity of sound. The formula is known to give qualitatively correct results 
when applied to materials under standard conditions,but its applicability can 
be expected to increase for materials subjected to high pressure. 

Finally,a few words about one more "influential" factor. It could be ques- 
tioned whether it is correct to compare equations which do not contain 
parameters of the same physical nature:equation (4) contains the temper- 
ature,while eq.(6) does not. It is true that the SK theory has been developed 
for the case T=OK.However,this should be understood just as a simplify- 
ing assumption, whose physical meaning is that this theory is applicable to 
materials at small teemperatures and subjected to high pressure,for which 
ep/kBT >- 1 (see Eliezer et al. ,1986 for details). 



4. Conclusions 

In this letter we have discussed a simple method for determining the 
thermal component of the EOS of solids under high pressure, thus correcting 
a small error made in a similar discussion (Celebonovic, 1982). Application 
to the interior of the Earth gives results in good agreement with geophys- 
ical data. Possible influence of several factors on the values obtained has 
been described in some detail. Using this method in laboratory high pressure 
work would necessitate an improvement in the calculation of the velocity of 
sound, and the intorduction of a suitable density-dependent cut-off frequency 
in equation (3). 
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